PARABOLIC EQUATIONS WITH DYNAMICAL BOUNDARY 
CONDITIONS AND SOURCE TERMS ON INTERFACES 



A.F.M. TER Elst, M. Meyries, and J. Rehberg 

Abstract. We consider parabolic equations with mixed boundary conditions and do- 
main inhomogeneities supported on a lower dimensional hypersurface, enforcing a jump 
in the conormal derivative. Only minimal regularity assumptions on the domain and the 
coefficients are imposed. It is shown that the corresponding linear operator enjoys max- 
imal parabolic regularity in a suitable L p -setting. The linear results suffice to treat also 
the corresponding nondegenerate quasilinear problems. 



1. Introduction 

In this article we are interested in the linear parabolic initial-boundary value problem 



ed t u - V • /iVu = / n in Jx(fi\S), (1.1) 

u = on Jx (dn\T), (1.2) 

ed t u + v ■ /xVtt + bu = fr on J x T, (1.3) 

ed t u + • /iV-u] = / s on J x E, (1.4) 

u(0) = m in fiUT, (1.5) 



and in its quasilinear variants. Here J = (0,T) is a bounded time interval, C is a 
bounded domain, T C <9f2 is a part of the boundary with outer normal u, and E C is 
e.g. a finite union of hyp ersurf aces, equipped with a normal field u-s. By ■ /xVtt] we 
denote the jump of z/s • fiVu over E. The case that T or E is an empty set is not excluded. 
We treat a nonsmooth geometry; e.g., it suffices that T and E satisfy certain Lipschitz 
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conditions. Nothing is supposed on the Dirichlet part dQ \ T of the boundary, and the 
boundary parts T and dQ \ T are allowed to meet. 

Also on the coefficients we impose only low regularity conditions. The (possibly non- 
symmetric) coefficient matrix /i is bounded and uniformly elliptic, e is positive, bounded 
and bounded away from zero, and b only has to live in an L p -space. The (possibly nonau- 
tonomous) inhomogeneities fa, fr, /s and the initial value w are assumed to be given. 

Parabolic problems with dynamical boundary conditions are considered by many au- 
thors, see e.g. |AmE] . |Escj . |Hin] |AQRB| , |BBR] and |BC] . but there always severe as- 
sumptions on the data, as smoothness, are imposed (compare also |FGGR] and |VVj , where 
the boundary condition on J x T is understood as Wentzell's boundary condition). It is 
the aim of this work to show that any smoothness assumption on the domain and the 
coefficient function // can be avoided. In particular, the domain Q does not need to be 
a Lipschitz domain. Let us briefly comment on this: a moment's thought shows that, by 
far, many natural domains fail to be Lipschitzian. For example, if one removes from a ball 
one half of its equatorial plane, then the remainder fails to be Lipschitzian. As another 
example, consider a pair of pincers as in Figure 1. It is also not Lipschitzian. The crucial 



Figure 1. A pair of pincers is not a Lipschitz domain. 



point is that such objects, obviously, occur in the physical world. In this paper we also 
allow the inhomogeneities not only to live in the volume of the domain, but to incorporate 
a part which is supported on the set E of lower dimension d — 1. This largely extends 
the applicability of the theory to real- world problems. The reader may think, e.g., of a 
heat source which is concentrated on an interface. Alternatively, one meets such constel- 
lations in electricity: surface charge densities induce a jump in the normal component of 
the dielectric displacement, see e.g. |Tam| Chapter 1]. 

Our approach to (1.1 )— (1.5 ), which also covers the case that T or S is empty, is essentially 
based on the theory of sesquilinear forms and the suitable incorporation of the boundary 
conditions into an LP- space. 
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We consider the approach in more detail. The boundary part T is Lipschitz regular, and 
the interface £ C £1 is a (d — l)-set in the sense of Jonsson-Wallin |JW] (cf. Assumptions 



2.2 and 2.4). For the equations we first treat the case e = 1, and consider the sesquilinear 



form 



t[u, v] — / [iVu-Vvdx, 
Jn 



which is defined on the space W r ' of W ' (f2)-functions vanishing on d£l \ T. Note that 
this reflects the Dirichlet conditions. For all u G W^ 2 we define the trace trtt on Y U £ 
in a suitable sense (based on |JWj ). and show that the map $u = (u,tru) is continuous 
and has dense ranee from Wl' 1 into L 2 := L 2 (VL) ® L 2 (T U E; dU d -i) (see Lemma |2.10[). 



Here Hd-i is the (d — l)-dimensional Hausdorff measure. These properties of the trace are 
a consequence of the regularity of V and E. As the form t satisfies an ellipticity condition 
with respect to Z, the results in |AEj imply that t induces an operator A 2 on L 2 . For all 
<f G W^ 2 and $<p G Dou^^) its constitutive relation is 

/ {A 2 Z<p)3${dx + dU d -x) = f nV<p- V^dx, f6^ 2 ((l). (1.6) 
Jnurws Jn 



Let us show that A 2 describes the spatial derivatives occurring in (1.1), (1.3) and (1.4), 
respectively, in an adequate manner. The argument is heuristic in general; moreover we 
identify within these calculations (p with 3y2, in order to make the writing more suggestive. 
Let A be a surface which is piecewise C l and which decomposes Q into two subdomains 
Qi and Q 2 . (A prototypical situation is when Q is a circular cylinder, T is its upper plate, 
and E is the midplane of fL) First put E = A n and assume that the outer normal v\ of 



Qi across A equals v% on dfli n S. According to (1.6), for all cp G Dom(v4 2 ) we have 



/ (A 2 <p) ip (dx + dHd-i) = j fiV(f ■ Vip dx + /iVy? • dx 
inurus Jn x Jn 2 

for all ip e W^ Z (Q). Since ip vanishes on dfl \ T, one can apply Gauss' theorem to obtain 

/ fiVip-Vtp dx = / (—V-fiVip)ipdx+ (v-fiVip) ip dHd-i+ / (z/r/iVy?) ip dH d -\- 
Jn x Jn x Jdn x v\Y J Ann 

(1.7) 



An equation, analogous to (1.7), can also be written for f2 2 . Then the unit normal v 2 of 



Q 2 across A equals —v\ and one deduces 

/ (A 2V )^ (dx + dHd-i) = / (-V ■/iV<^)V'Gte+ (u- fiV^dUd-i (1.8) 
inurus Jn Jv 

[vy, ■ nVflipdWd-i, 

Ann 
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where [z/ s • /xVy] = z/ s • (/iV<^|on in s — A^^l^ans) is the jump in the conormal derivative. 



Thus, varying ip suitably and comparing both sides of (1.8), one recognizes that A2 has in 
fact three 'components', namely 

(1) the divergence of the vector field /xV<£> on Q \ E, taking L 2 (f2)-functions as values; 

(2) the conormal derivative on T, taking L 2 (T; d7^d_i)-functions as values; 

(3) the jump in the conormal derivative on E, taking L 2 (E; tfHd_i)-functions as values. 



If one takes E as a proper subset of A fl (which admits the (d — l)-property), then (1.8) 
leads to the equation 



/ {A 2 ip) ip dU d -i = [z/ s ■ /iV(p] if) dHd-i, 
Jt, JAnn 



which enforces [u% ■ /xVy] to vanish on (An Q) \ E. Hence the dynamic equations on V and 
E are modelled by the part L 2 (ruE; dHd-i) of the base space L 2 . The subsequent analysis 
will show that, in either the elliptic or in the parabolic setting, these three components may 
be prescribed, and the equation indeed has a solution in the functional analytic setting 
which we will establish. Moreover, the solution depends continuously on the data. 

The operator — A2 generates a holomorphic, submarkovian Co-semigroup of contractions 
on L 2 , and may thus be extended to a semigroup of contractions on L p for all p6 [1, 00]. 
Denoting the corresponding generators by — A p , it turns out that for all p G (1, 00) the oper- 
ator —e^ l A p generates a holomorphic Co-semigroup of contractions on a suitably renormed 
L p -space. This has two important consequences. First, applying an abstract result that 
is presented e.g. in [LXL Proposition 2.2], we obtain a bounded holomorphic functional 
calculus for e~ 1 A p with angle strictly smaller than |, and in particular the boundedness 



of the purely imaginary powers (see Theorem 2.22). Moreover, the pioneering theorem of 
Lamberton |Lamj gives us maximal parabolic regularity for e~ l A p in Theorem 3.4, which 
we consider as the main result of this work. The introduction of temporal weights as in 
jPSj further allows to reduce the regularity of the initial data almost up to the base space 



L p . This yields the solution of (1.1 )— ( 1.5 ) in an adequate manner, see Theorem 3.10 

Based on these linear results we treat a nondegenerate quasilinear variant of (1.1 )— (1.5), 
even if the right hand side explicitly and discontinuously depends on time (Theorem 4.5). 
Here a difficulty is that the domain of the realization of the operator —V • /xV on L p is 
not independent of the coefficients \i. We therefore consider a problem which is obtained 
when applying the Kirchhoff transform to the original one, and which involves only one 



fixed operator (see Definition 4.2). Maximal parabolic regularity then allows to apply a 
result of Priiss |Priij (see also |C1LJ ) to the transformed problem, giving local existence and 
uniqueness of solutions in a suitable sense. Throughout it is essential that Dom(Ap) C L°° 
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for large p and 6 sufficiently close to f , which is a consequence of ultracontractivity esti- 



mates for the semigroup (see Lemma 2.19 and Proposition 3.5). The quasilinear problems 
may be of relevance for the applications: the heat source on the hypersurface can depend 
on the solution itself, and, additionally, explicitly on time. 

Let us briefly compare the approach in this paper with those in [Gri2j . [HaRlj and 
|HaR3] for static Robin boundary conditions. There the Banach space under consideration 
is a negatively indexed Sobolev space of type H~ e,q or a Sobolev-Morrey space. In contrast 
to that settings, in L p one may form the dual pairing of the above parabolic equation with 
the indicator function xa of suitable subsets A C Q. Then one may, additionally, apply 
Gauss' theorem to (— V ■ fiVu, Xa) — J A — V • fiVu dx + J AnS —V ■ /iVw d%d-i- This allows 
to recover the underlying physical balance law for the parabolic equation, which is the 
starting point for the numerical treatment of such problems. For more details we refer to 
Remark 13.131 

This paper is organized as follows. In Section [2] we introduce the spaces L p , define 
an appropriate realization of —V • /iV and show that it admits a bounded holomorphic 



functional calculus. In Section [3] we show that in this setting (1.1 )— ( 1.5 ) enjoys maximal 
parabolic regularity, and in Section [4] we treat the quasilinear case. We finish with some 
concluding remarks in Section [5] 



2. Elliptic operators on L p 

2.1. Notation. Throughout this paper C(X;Y) denotes the space of bounded linear op- 
erators from X to Y, where X and Y are Banach spaces. If X — Y, then we abbreviate 
C(X). Note that if X and Y are two Banach spaces spaces such that Ic^as vector 
spaces, and both X and Y are continuously embedded in a Hausdorff locally convex space, 
then the inclusion map from X into Y is continuous by the closed graph theorem. 

In the sequel let Q be a bounded domain in Mr with d > 1 and F an open part of its 
boundary dfl, which may be empty. If p 6 [l,oo), then L P (Q) is the space of complex- 
valued, Lebesgue measurable, p-integrable functions on Q, and for all 6 G [0, 1] we denote 
by W e,p (Q) the usual Sobolev-Slobodetskii spaces, see |Gri] or |Mazj . Moreover, L°°(Q) 
is the space of Lebesgue measurable, essentially bounded functions on Q. The (d — 1)- 
dimensional Hausdorff measure on M. d is denoted by %d-i- We denote by B(x,r) the ball 
in ]R d centred at x with radius r. 



2.2. The function spaces. In this subsection we consider the function spaces on which 
(1.1 )— ( 1.5) will be posed. 
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Definition 2.1. For all q G [1, oo] we define W^ 9 as the closure in W l,q (VL) of the set 

c?°(ft) d = : u e c™{R d ), supp(m) n {an \ r) = 0}. 

Throughout this paper we make the following assumption on T. 

Assumption 2.2. For all x G V there is an open neighbourhood V x of x and a bi-Lipschitz 
mapping F x from V x onto the open unit cube E in M. d , such that -F x (:r) = and F X (Q fl V x ) 
is equal to the lower open half cube E_ = (— 1, x (—1, 0) of E. 

The reader should notice that the domain Q does not need to be Lipschitzian. Moreover, 
nothing is supposed on the boundary of T within dfl. 

An important technical tool is an extension operator for the Wj^-spaces. 

Proposition 2.3. There is an extension operator (£: — > L 1 (IR d ) such that the re- 
striction <E\ w i,q maps W^ 9 continuously into W 1,q (M. d ) for all q G [l,oo], the restriction 
rnaps L q (Q) continuously into L q (M. d ) for all q G [l,oo] and supp <£u C B(0,2R) 
for all u G L 1 ^), where R = sup{|x| : x G fl}. 

Proof. The proof is given in |ER4 Lemma 3.4] for the case q = 2, but carries over to all 
q G [l,oo]. Moreover, the second assertion is also easily checked. The last statement 
follows by multiplication with a suitable (M d )-function. □ 

It turns out that a classical condition from geometric measure theory is tailor made in 
order to define a geometric assumption on a (d — l)-dimensional shape S in Q. 

Assumption 2.4. Let X C Q be a (d — l)-set in the sense of Jonsson-Wallin \JW\ Sub- 
section VII. 1.1]. Precisely: the set S is Borel measurable and there exist Ci,c 2 > such 
that 

cir 4 - 1 < n d -i(B(x,r) n S) < c 2 r d - 1 (2.1) 
for all x G S and r G (0, 1). 

Remark 2.5. We emphasize that S does not have to be closed. Nevertheless has £ 
finite (d — l)-dimensional Hausdorff measure, according to (2.1). The prototype of S is 
the finite union [J . Sj of Lipschitzian hypersurfaces. In that case the restriction of the 
Hausdorff measure Hd-i to T or to Sj can be constructed explicitly in terms of the local 
bi-Lipschitz charts (compare |EGl Section 3.3.4 C]). In particular, if S is a finite union of 
Lipschitz graphs, then (2.1) is easily verified using this representation of TLd-i- Moreover, 
Assumption 2.4 implies for general E that E is of (c/-dimensional) Lebesgue measure 0. 



Throughout this paper we always presume Assumptions |2.2| and 2.4 
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Definition 2.6. We denote by p the restriction of the Hausdorff measure Hd-i to TU E. 
If u G L 1 1 oc (M d ) and F C W 1 is a set, then define the function tipu as in |JW| Page 15] 

by 



(tr F u) (x) = lim — / u(y) dy, 

r->o \B(x,r)\ J B{XiT) 

for all x G F for which the limit exists. The domain Dom(trp-a) of txpu is the set of all 
x G F for which this limit exists. 

Lemma 2.7. Let q,r G [1, oo) and 9 G [0,1]. Let (£ be the extension operator as in 



Proposition 2.3 



(i) If--^f-<l, then <E maps W^' q continuously into W e ' r (R d ). 

(ii) If - — < i, i/ien <£ maps Wp' 9 compactly into W 9 ' r (R d ). 

(iii) 7/0 G (-, 1], i/ien ine trace map u i— >■ trrus"a is continuous from W e,q (W 1 ) into 
L"(ruE;(ip). 



Proof. 1 (i) 1 and ' (ii) . This follows from Proposition 2.3 , the support property of <E and the 
usual Sobolev embedding. 

Since T and £ are disjoint, the natural map from the space L q (T U S; dp) into 



m 



L q (H; dHd-i) x L q (T; dHd-i) is a linear, topological isomorphism. Therefore, it suffices to 
show that the trace maps u h-> trrw and -a i— > tr^-a are continuous from W ,q (M. d ) into 
L q {T;dU d -x) and a"H d _i). 



It follows from [JW| Chapter VIII, Proposition 1] that property (2.1) inherits to the 
closure £ of £. Then the trace operator u \- > tr^u is bounded from W e ' q (W i ) into 
L 9 (S; dHa-i) by |JW[ Chapter V, Theorem 1]. But the set difference S\E is of Hd-i mea- 
sure (see again |JWl Chapter VIII, Proposition 1]). Consequently the spaces L q (T>; d%d-i) 
and L q (Y>; dUd-i) are identical. 



Next we consider the set T. Using the notation as in Assumption |2.2| , for every x G T 
the map F x provides a bi-Lipschitz parametrization of dfl (1 V x , where the parameters 
run through the upper plate P := (—1,1) x {0} of the half cube EL. Moreover, the 
Hausdorff measure H d -\ on dQ fl V x is the surface measure, and the latter is obtained 
from the Lebesgue measure on (—1, l) ^ 1 x {0} via the bi-Lipschitzian parametrization, 
see (ESI Section 3.3.4 C]. Define W x = F x ((-\, l)*" 1 x {0}). Then W x C dQ. There 
exist n G N and X\, . . . ,x n G T such that W Xl , • • • , Wx n is a finite cover of T. Obviously, 
W Xl , • • • , yVx n is also a finite cover of T. Moreover, it is not hard to see that Uj=i W Xj 
is a (d — l)-set in the sense of Jonsson-Wallin (compare |HaR2l Lemma 3.2]). Hence 
by |JW| Chapter V, Theorem 1] there exists a continuous trace operator from W ,q (M. ) 
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into L 9 (U" =1 W a:j ; dHd-i)- Combining this operator with the restriction operator to T, one 
obtains the desired trace operator into L q (T; dWd-i)- D 

For all u G L\ oc {VL) define the function tiu as in |JWl Section VIII. 1.1] by 

Dom(tr M ) = L e r U £ : lim ( u(y) d y ex,st 8 } 

and ^ 

(trw)(x) = lim — — — / u(y) dy 

r^o \B(x,r) niL\ J B { x , r )m 

for all x G Dom(tru). 

The above defined trace enjoys the following mapping properties. 

Proposition 2.8. Letq,r G (1, oo) and suppose that ^ < Thentiu G L r (TUT,;dp) 
for all u G W^' q , and the map u i— >• tr it zs compact from W^ q into L r (T U E; dp). 



Proof. Let £ be the extension operator as in Proposition |2.3[ Then it follows from 
Lemma 2.7 that u i-)- trrus maps Wp' 9 compactly into U\Y U E; dp). But if u G Wr' 9 , 



then we claim that 

(trM)(x) = (trru S Cu)(ar) (2.2) 
for "Hd-i-a.e. x G ruE. Obviously, this identity holds for "H^i-a.e. x G E since Scfi. For 
"Hd-i-a.e. x G T we can argue as in the proof of |JW| Chapter VIII, Proposition 2], where 
the case V = dQ is considered. Indeed, the arguments given there are purely local. Since 
<Eu G W 1 '^^) it follows that for U d G T there exists a Borel set £7 C M. d such 

that ^_i(£nB(i,r)) = o^r^ 1 ) and (<£u)(x) = lim (<Eu)(y). Using these properties 

y-rx, y$E 



of E, the same arguments as in the last part of the proof given in |JWj establish (2.2). □ 



The space on which (1.1 )— ( 1.5 ) will be realized is given as follows. 

Definition 2.9. For all p G [1, oo], denote by L p the Lebesgue space L p (fl U T; dx + dp). 
We denote the space of all real valued functions in L p by L^. 

Observe that there is a natural topological isomorphism between L p and the direct sum 
L P (Q) © L P (T U E; dp) and we will identify W with L P (Q) © L P (T U E; dp) through this 
natural map. 



By Proposition 2.8 we can define the map W^' 2 — > L 2 by 

Zu = {u, tr u) G L 2 (tt) © L 2 (T U E; dp) = L 2 . 



Note that one can always choose some p > 2 in Statement (ii) of the next lemma. 
Lemma 2.10. 
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(i) The map 3 is continuous and has dense range. 

(ii) Ifp G [1, oo) and (d - 2)p < 2(d - 1) ; then 3Wp 2 C U 

(iii) The map $ is compact. 



Proof. 
2 



i)'. The continuity follows from Proposition 2.8 Let / 



(/n, /a) e L 2 (Q)® 

L 2 (T U E; dp) and suppose that (3m, /)i 2 (n)eL 2 (rus ; dp) = for all u G Wjl' 2 . We show that 
/ = 0. For all u G C c °°(fi\E) one has = (5m, /) = J Q ufadx. Since L7 c °°(fi\E) is dense in 
L 2 (fi \ S) = L 2 (Q) one deduces that fa = 0. Therefore = J ruE trM /a dp for all w G H^p' 2 
and in particular for all u G Cp°(f2). But {w|ru£ : w G C r <> (r2)} is dense in L 2 (T U E; dp). 
So / 9 = 0. 

' (ii)[ . If (£ is the extension operator as in Proposition 2.3| then it follows from Lemma 2.7|[i) 
that (£ maps W^' 2 continuously into L p (IR d ) for all p G [l,oo) with (d — 2)p < 2d. So 



W^ 2 C L P (Q). Now the statement follows from Proposition 



m 



2.1 



.... 



2.8 



that the restriction (B\a maps 



1.2 



It follows immediately from Lemma 
compactly into L 2 (Q). So the embedding of W^' 2 into L 2 (Q) is compact. Also the map 
tr is compact from W^ 2 into L 2 (T U E; dp) by Proposition 
compact. 

We end this subsection with a truncation lemma. 



2.8 



Therefore the map Z is 

□ 



Lemma 2.11. Let u G W r ' be real-valued. Then u A In G VF r ' and 3(w A In) 
(5m) A Inur- 



Proof. The first statement is shown in the proof of |ERl Theorem 3.1]. The second 
statement is obvious for real- valued u G C^(Q). Since the maps u H- Z(u A In) and 
u i — y ($u) A Inur are continuous on the real version of W^' 2 , the identity carries over to the 
general case by density. □ 

2.3. The operator on L p . In this subsection we introduce a differential operator on L p 



that corresponds to the spatial derivatives in (1.1), (1.3) and (1.4). 



Throughout the remaining of this paper we adopt the next assumption. 

Assumption 2.12. Let p = {p^} : ^ — >■ £(1^1^) be a measurable map from fl into 
the set of real d x d matrices. We assume that there are p,,p* > such that 

d d 

\\K x )\\cW d ) < and Vk,i( x )£k£i > At. 5Z^fc 



k,l=l 



k=l 



for all x G f2 and £ = (£i, . . . , £d) G 



We emphasize that p does not have to be symmetric. 
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Definition 2.13. Define the sesquilinear form t: W r ' x W r ' — > C by 

t[u, v] — / [iVu-Vvdx. 
Jq 

1 2 

We emphasize that the domain of the form t is the space W r ' , which appropriately incor- 
porates the Dirichlet condition on <9f2\r, compare [CiaJ Section 1.2] or |GGZl Section II. 2]. 
The form t is continuous and 

Ret[u,u] + \\M\h > (A*- a !)IMIw l . a ( 2 - 3 ) 



for all u G W r ' . Therefore by Lemma 2.1C i) and |AEl Theorem 2.1] there exists a unique 



operator A 2 in L 2 such that for all ip,ip G L 2 one has if G Dom(A 2 ) and A 2 ip = ip if and 

12 

only if there exists au£ W r ' such that Zu = tp and 

t[u,v] = (ip,3v) h 2 (2.4) 

for all v G W^ 2 ■ Although the form domain of t is W^' 2 , the operator A 2 is an operator 
in L 2 . We refer to the introduction for a discussion of the relation of A 2 to the original 
problem (1.1 )— ( 1.5 ). 

Remark 2.14. The construction of A 2 generalizes the derivation of an operator from a 
suitable form s to the case when the form domain D s is a priori not contained in the 
corresponding Hilbert space (compare |Katl Section VI. 2] for the classical case). The 
substitute for the inclusion D s C is the definition of an appropriate embedding operator 
Z'- D s — > $). Fortunately, all tools for form methods are still available. 

Proposition 2.15. The operator A 2 is m-sectorial with vertex and semi-angle arctan 
Moreover, A 2 has compact resolvent. 

Proof. It follows from \AE\ Theorem 2.1] that A 2 is m-sectorial. Let ip G Dom(v4 2 ) and 
u G Wf 2 with %u = tp. Then Re(A 2 ip, </?) L 2 = Re t[u, u] > 0. Hence the vertex is 0. 
Further, one has Ret[w, u] > /i, J n \ V Re u\ 2 + \ V lmu\ 2 dx and 

|Imt[w,w]| < 2/i* / |VReu||VImu| dx < /i' / | V Re u\ 2 + | V Imw| 2 dx. 
Jq Jq 

Thus | a.rg(A 2 ip, (p)u | < arctan ^ if tp ^ 0. 



Since the map $ is compact by Lemma 2.1C 'hi) , the generator has compact resolvent by 



[AH Lemma 2.7]. □ 



We continue with the analysis of the operator A 2 . By Proposition 2.15 the operator A 2 
is m-sectorial with vertex and semi-angle arctan—. Hence by |Katj Theorem IX.1.24] 
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the operator —A 2 generates a holomorphic semigroup, denoted by S, which is holomorphic 
and contractive on the sector with semi-angle arctan — . 

Proposition 2.16. The semigroup S leaves invariant, it is submarkovian and positive. 

Proof. Clearly the set is closed and convex in L 2 . Moreover, (p i— y Re p is the projection 
from L 2 onto and Ret(u, u — Reu) = for all u G W^' 2 . Since the form t is accretive, 
the set Lj| is invariant under the semigroup by \AE\ Proposition 2.9(h)] 

Next, let C = {u G L 2 : u is real valued and u < 1}. Then C is closed and convex. Let 
P: L 2 — y C denote the orthogonal projection. Then Pu = (Rett) A Inur- Let u G Wjr. 



By Lemma 2.11 one has (Reu) A 1^ G W^' 2 and P$u = ^((Rett) A 1^). Moreover, an easy 
calculation gives 

Ret[(Rew) At n ,u- (Reu) A l n ] = 0. 
Observing that the form t is accretive, it follows from |AE[ Proposition 2.9(h)] that C is 
invariant under the semigroup S. Now let (p G L 2 D L°° and t > 0. There exists anaGl 
such that ||5t<^||L°° = || Re(e ia S'^)||L°°. But Re(e ia S t (p) = S t Re(e ia ip). Therefore 

HS^Hl- = \\S t Re(e"V)lk-> < || Re(e ta (p)\\ h oo < \\lp\\ l ™ 

and S is submarkovian. 

Finally, if p G and (p < 0, then mp G C for all n G N. So St(n<p) < 1 for all t > 
and n G N. Therefore St<p < and S is positive. □ 

Corollary 2.17. For all p G [1, oo] the semigroup S extends consistently to a contraction 
semigroup on h p . The semigroup is strongly continuous for all p G [1, oo). 



Proof. Observe that if the coefficient matrix jj, satisfies the conditions of Assumption 2.12 



then its transpose satisfies these as well. Thus the dual semigroup S* shares the same 
properties as S. Now the assertion follows from Proposition 2.16 and standard interpolation 
and duality arguments, see e.g. |Ouh2| page 56]. □ 

We denote the generator of by — A p . Then — A p is dissipative by the Lumer-Phillips 
theorem. If no confusion is possible we write S = S^'. 

Remark 2.18. It is possible to prove the dissipativity of — A p also by showing that the 
form — t is p-dissipative, cf. |CMj . 

Lemma 2.19. 

(i) The semigroup S is ultracontractive. Moreover, for all [3 > d — 1 and u > there 
exists a c > such that 

\\s t <p\\ hq < cr^H^ii^iiLP 
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fiv) 



for all t > , if G L p and p,q G [1, oo] with p < q. 

//l<P<g<oo and j G N are stzca i/iai (i — i) < 1 ; £/ien i/ie operator 

(A p + l) - -' maps L p continuously into h q . 

The operator A p has compact resolvent for all p G (1, oo). 

If the matrix of coefficients /i is symmetric, then the operator A 2 is self-adjoint and 
positive. 



Proof. Let r G (2, oo) be such that (d — 2)r < 2(d — 1). Then it follows from 

t In i r\V--\\ ^1 X OTI/1,2 tt 7- 1 ,1 1 , 1 1 J 1 



Lemma 



2.1C 



that ZW T ' C L r , and the inclusion is continuous by the closed graph 



theorem. Let cp G L and t > 0. Since S^y? G Dorr^Aj), there is a u G W r ' such that 
Stf —Zu. For given w > one has 

\\SM\lr = \\M\lr <C\\u\\ 2 w i, 2 < C{n, A l) _1 (Ret[u, u] + \\3u\\h) 

= C(jjl. A l)- x ( Re(A 2 £^, S^) L 2 + ||%>||y 

<CrV<iM|£ 2 , 



for suitable C, C > 0, using (2.3), the definition of A 2 , the Cauchy-Schwarz inequality and 



the holomorphy and contractivity of St- Therefore the semigroup S is ultracontractive, and 
by |Ouh2t Lemma 6.1] there exists a c > such that 



-LOt 



<SW>||l°° < ct \\ip\ 



for all t > and <p G L 2 . Now duality and interpolation give Statement 
follows from 



Statement 



and the well-known formula 

1 \! / 



U 1)! 



t^e^Stdt. 



Statement 



m 



of the lemma is easy to verify. 



is a consequence of Proposition 2.15 and interpolation. The last statement 

□ 



2.4. Multipliers acting on Lebesgue spaces. In order to solve (1.1 )— (1.5 ), we divide 



(Li]) (at first formally) by e. Obviously, one is then confronted with the necessity to 
investigate the functional analytic properties of operators of the type qA p , where c is 
a bounded strictly positive measurable function. Concerning the generator property of 
an analytic semigroup in a space L P (Q) this was carried out in [UKRj and concerning 
maximal parabolic regularity on L P (Q) in |HiRj . In the latter case the decisive instrument 
was the insight from |DO] that a suitable multiplicative perturbation does not destroy 
upper Gaussian estimates, which in turn imply maximal parabolic regularity on L P (Q). 
Unfortunately, we cannot apply this here, since our Lebesgue space does not only live 'on 
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the volume'. But a surprisingly simple trick allows us to overcome the problem in the 
present context. 

The next proposition is of independent interest. 

Proposition 2.20. Let (X, B, A) be a measure space and letr: X — > (0, oo) be a measurable 
function such that both r and r _1 are bounded. Let p G [1, oo) and let T be an operator in 
D>(X, dX). 

(i) If T is dissipative on L p (X,dX), then tT is dissipative on L p (X,r~ l dX). 

(ii) If T generates a strongly continuous contraction semigroup on L p (X,dX), then tT 
generates a strongly continuous contraction semigroup on L p (X,r~ l dX) . 

(iii) If p = 2, 6 £ (0, |] and T generates a holomorphic semigroup in L 2 (X,dX) which is 
contractive in the sector with semi-angle 9, then tT generates a holomorphic semi- 
group in L 2 (X,r~ 1 dX) which is contractive in the sector with semi-angle 9. 

Now suppose that p = 2 and T generates a strongly continuous contraction semigroup S 
on L 2 (X,dX). Denote the semigroup generated by tT on L 2 (X,r~ 1 dX) by S T . 

(iv) // S leaves the real valued functions invariant, then S T also leaves the real valued 
functions invariant. 

(v) If S is positive, then S T is also positive. 

(vi) Suppose S is submarkovian. Then S T is also submarkovian. Hence for all q G [2, oo) 
the semigroups S and S T extend consistently to a strongly continuous semigroup 
and S( T ' q ^ on L q (X,dX) and L q (X,r~ 1 dX), respectively. Let T q and T Tq denote the 
generators. Then T Tq = rT q for all q G [2, oo). 



Proof. |(i)[ . The operator T is dissipative on L P (X, dX) if and only if 

Re / (Tf)\f\ p - 2 JdX<0 

for all / G D(T). This implies the dissipativity of tT on L p (X,r~ 1 dX). 

(n)j. Since T generates a contraction semigroup on L p (X,dX), it follows that T is 



dissipative. Therefore tT is dissipative on L p (X,r 1 dX) by Statement (i) So by the 
Lumer-Phillips theorem it remains to show that the operator tT — 1 is surjective on 
L p (X,r' 1 dX). 

Let 5 > be such that r _1 — 5 > 0. Then the multiplication operator — (r _1 — 5) is 
dissipative on L P (X, dX) and has a relative bound equal to zero with respect to T. Therefore 
T — (t -1 — 5) generates a strongly continuous contraction semigroup on L P (X, dX) by the 
perturbation result |Davt Theorem 3.7]. Hence T — r" 1 is surjective on L p (X,dX) by the 
Lumer-Phillips theorem. But this implies that tT — 1 is surjective on L p (X,r~ 1 dX). 
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For all a G (— 9, 9) the above applies to the operator e ia T . Therefore e ia rT 
generates a strongly continuous contraction semigroup on L 2 (X,r~ 1 dX). Hence by \K&t\ 
Theorem IX.1.23] the operator tT generates a holomorphic semigroup in L 2 (X, r~ 1 d\) 
which is contractive on the sector with semi-angle 9. 

Now suppose p = 2 and T generates a strongly continuous contraction semigroup S on 
L 2 (X, dX). Let C be a closed convex subset of L 2 (X, dX). Then C is also closed and convex 
in L 2 (X, r~ l dX). Since T is m-dissipative it follows from |Ouhl[ Theorem 2.2] that C is 
invariant under S if and only if Re(T/, / — Pf)i, 2 (x,d\) < for all / G D(T), where P 
is the orthogonal projection in L 2 (X,dX) onto C. Similarly, since tT is m-dissipative on 
L 2 (X, r~ 1 dX), the set C is invariant under S T if and only if Re(rT/, f — P T f) L 2 {x,r- 1 d\) < 
for all / G D(tT), where P T is the orthogonal projection in L 2 (X,r~ 1 dX) onto C. But 
D(tT) = D(T). Hence if P = P T , then C is invariant under S if and only if C is invariant 
under S T . 



Then for the proof of Statement (iv) choose C = {/ G L 2 (X, <iA) : / is real valued} 



and note that the projection is Pf = Re/ = P T f. For the proof of Statement 
choose C = {/ G L 2 (X, c?A) : / is real valued and / > 0} and note that the projection is 



Pf = (Re/) + = P T f. For the submarkovian part in the proof of Statement (vi) choose 
C = {/ G L 2 (X,dX) : \ f\ <l a.e.} and note that the projection is Pf = (\f\ A l)sgn/ = 
P T f- 



It remains to prove the second part of Statement (vi) Let q G [2, oo). Let u G D(T Tq ) fl 
D{T T}2 ). Write v = T Ty2 u. Then u G L 2 (X,dX) n L q {X,dX) and T Tj9 m = T Ti2 m = u. So 
u G L q (X,dX) and r -1 !! G L q (X,dX) since r -1 is bounded. Moreover, T r> 2 = TT2, so 
u G DiT?) and T 2 m = r _1 f . Therefore 



t-\S^u -u) = t~\Sl 2) u -u) = r 1 [ Si 2) T 2 uds = r 1 [ Si q ~>T 2 uds 

Jo Jo 

for t > 0. As £ 4- 0, the latter term converges to T 2 u in L q (X, dX) by the strong continuity 
of S'W. Hence it G D(T q ) and T g w = T 2 n = r _1 t>. Then tT^m = t> = T r (J M. We proved that 

D(T T>q ) n D(T r , 2 ) C D(rT g ) n ^(rT,) 

and T T (? u = rT g M for all u G D(T Tjq ) fl D(T Tj2 ). Similarly the converse inclusion is valid, so 

D{T TS ) H D(r T)2 ) = D{rT q ) n ^(rT,) = D(T,) fl D(T 2 ). 

We claim that D{T q ) fl D(T 2 ) is dense in D(T q ) = D{rT q ). Consider the set 

X? = {r 1 / S-^uds : it G L q {X,dX) nL 2 {X,dX),t G (0,oo)}. 
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Then V C D(T q ). Since and are consistent, also V C D(T 2 ). So V C D(T ff ) n 
D(T 2 ). Moreover, lim^o*" 1 jj Scurfs = it in L 9 (X, dA) for all u G L 9 (X, dA) n L 2 (X, dA) 
and L«(X, dA)fU 2 (X, dA) is dense in L q (X, dX). Therefore V is dense in L q (X, dX). Clearly 
T> is invariant under S^ q \ Hence T> is a core for T 9 by |EN] Proposition II. 1.7]. This implies 
that D{T q ) fl D(T 2 ) is dense in D{T q ). The same arguments show that D(T T}<] ) fl D{T T ^) is 
dense in D(T Ttq ). Hence T Tj(? = rT g . □ 

Let <r: f2 U T — 7- (0, oo) be a measurable function such that q, q^ 1 G L° 



We write 



:= LP(fiur;? _1 (da: + dp)). 



Proposition 2.20 allows to transfer the conclusion of Corollary 2.17 to the operators qA p . 

Theorem 2.21. For all p G [1, oo) the operator — qA p generates a strongly continuous 
positive semigroup S^ ,p} of contractions on the space L|. The semigroups are consistent. 
Moreover, S^' 2 ^ is holomorphic and contractive in the sector with semi-angle arctan—. 



Proof. For p > 2 all follows from Propositions |2. 15] |2 . 1 6| and |2 .20| The dual of the operator 
qA 2 on L 2 is given by <rA#, where A* is the operator obtained with coefficient matrix equal 

the dual semigroup (S 1 * 5,2 ))* 



2.20 



to the transpose of the matrix /x. Hence by Proposition 
is submarkovian and extends consistently to a strongly continuous semigroup on for 
all p G [2,oo). Then by duality the semigroup S^' 2 ^ extends consistently to a strongly 
continuous semigroup on L£ for all p G [1, 2]. □ 

2.5. Consequences for the operators qA p on L p . We have the following abstract prop- 
erties for <;A p . 

Theorem 2.22. Let p G (1, oo). Then one has the following. 

(i) The operator qA p admits a bounded holomorphic functional calculus on h p , with angle 
strictly smaller than | . In particular, it admits bounded imaginary powers. 

(ii) For all 9 G (0, 1) one has 



71 



HI 



If 6 G (0,1], then Bom((qA p ) e ) = [L p , Dom(qA p )] 
complex interpolation. 



poo 

/ t-\qA p + 1 + ty 1 dt. 
Jo 

Dom(Ap), where [•, -]g denotes 



Proof, (i) 1 . For all p G [1, oo) denote by S^'^ the contraction semigroup on L£ generated 
by — qA p . Then the semigroups 

S^ri w ith P e [1 , oo) are consistent. Moreover, S^' 2 ^ 
is holomorphic and bounded on a sector. Let p G (l,oo). Then it follows from |Ouh2t 
Proposition 3.12] and duality that S^' p ^ is holomorphic and bounded in a sector (which 



16 A.F.M. TER Elst, M. Meyries, and J. Rehberg 

depends on p). Also S^'^ is a positive contraction semigroup. Hence the operator S^'^ is 
contractively regular for all t > 0. So by |LX[ Proposition 2.2] the operator <;A P admits a 
bounded holomorphic functional calculus on L£, with angle strictly smaller than |. This 
is then also the case on L p , since h p = as vector spaces, with equivalent norms. 

(ii) [ . For the integral representation see |Lun[ (4.41)]. 

(iii) . Since qA p admits bounded imaginary powers, it follows from |Lunl Theorem 4.17] 



that 

Dom((^y e ) = [L p ,Dom(^4 p )] e . 

Since Dom(?i p ) = Dom(A p ), one has Dom((sA p ) e ) = [L p , Dom.{A p )\e, and the result 
follows. □ 

3. Linear parabolic equations 
In this section we will draw conclusions for linear parabolic equations, which, in partic- 



ular, allow to give (1.1 )— ( 1.5 ) a precise meaning and afterwards to solve it. 

In the following, J = (0, T) denotes a bounded interval and X a Banach space. Through- 
out we fix the numbers 

1 < s < oo and - < a < 1. 

s 

We introduce the weighted space 

L s a (J;X) = {u: J X : [t i-> t l ~ a u{t)} G L S (J;X)}, 

and the corresponding weighted Sobolev space 

W l a 3 (J-X) = {u e L s a (J-X) : v! e L s a (J;X)}, 

where here and below the time derivative is understood in the sense of X-valued distribu- 
tions (see |Ama| Subsection III. 1.1]). These are Banach spaces when equipped with their 
canonical norm, respectively. Note that a = 1 corresponds to the unweighted case, i.e., 
L{ = L s . By [FSl Lemma 2.1] one has W^> S (J;X) C W l ^(J;X), which implies that each 
element of W^ ,S (J; X) has a well-defined trace at t — 0. 

Definition 3.1. Let A be a closed linear operator on X with dense domain Dom(A). We 
say that A has maximal parabolic L s a (J; X) -regularity, if for all / G L s a (J;X) there is a 
unique solution u G W^' S (J; X) n L s a (J; Dom(v4)) of 

u + Au = f, u(0) = 0. 

We write MR s a (J]X) for the class of all operators on X with this property. 
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We proceed with some comments concerning maximal parabolic regularity. 

(1) It is shown in [PHI Theorem 2.4] that A E MR{{J\ X) if and only if A E MR s a (J; X) 
for all a E (-, 1], i.e., maximal parabolic L*-regularity is independent of the weight. 
(In fact, in [PS] only the case J = (0, oo) is treated, but the arguments given there 
also apply to bounded J.) In this sense it is natural to consider the temporal 
weights in the context of parabolic problems. 

(2) If A E MR{°(J - X) for an interval J = (0, T ), where T E (0, oo) and s E (1, oo), 
then A E MR s a (J; X) for all T <G (0, oo), s E (1, oo) and a E (J, 1]. This is shown 
in |Dorl Corollary 5.4 and Theorem 7.1]. In this spirit, we then simply say that A 
satisfies maximal parabolic regularity on X. 

(3) The notion 'maximal parabolic regularity' does not depend on the concrete norm 
of the Banach space. In other words: an operator A, satisfying maximal parabolic 
regularity on X, remains to satisfy maximal parabolic regularity if X is equipped 
with an equivalent norm. 

(4) If A satisfies maximal parabolic regularity on X, then —A generates an analytic 
Co-semigroup (cf. |Dorl Corollary 4.4]). If X is a Hilbert space, then the converse 
is also true, cf. [DcSJ. 

For the case of nontrivial initial values, the following has been proved in \PS\ Theo- 
rem 3.2]. We denote by (-,-)e,s the real interpolation functor, cf. |Tril Sections 1.3 and 
1.6]. 

Proposition 3.2. Suppose that A satisfies maximal parabolic regularity on X . Then for 
all f E L s a (J; X) and u$ E (X, Dom(v4)) Q _i the Cauchy problem 

s ' 

u + Au = f, u(0) = uo, 
has a unique solution u E W^' S (J; X) D L s a (J; Dom(A)), and the estimate 

\\u'\\l^(J;X) + ||w|U»(J;Dom(i4.)) < c(||uo||(X,Dom(A)) _i + (3-1) 

is valid for some constant c, independent of f and uq. 

By working in temporally weighted spaces one can thus reduce the regularity of the 
initial values u$ almost up to the base space X. 

We have the following embeddings for the weighted maximal regularity class. The space 
of 7-Holder continuous functions is denoted by C 1 . 

Proposition 3.3. If A satisfies maximal parabolic regularity on X, then 
W l a s {J- X) n L s a (J; Dom(A)) C BUC(J; (X, Dom(A)) Q _i J n C(J; (X, Dom^))^! J. 

s 5 s ' 
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Moreover, for every 9 G [0, a - i) there is a 7 G (0, 1) such that 

W^ s (J;X)nL s a (J;Dom(A)) C C"V; [X,Dom(A)) ). 

Proof. The first inclusion is shown in |PS| Proposition 3.1]. The second one can be proved 
along the lines of |DMRTt Lemma 1]. □ 

We apply a classical result of Lamberton |Lamj to the operators qA p . 

Theorem 3.4. Let q: Q U T — > (0, oo) be a measurable function such that G L°°. 

Then for all p G (1, oo) the operator qA p satisfies maximal parabolic regularity on L p . 



Proof. Theorem 2.21 states that the semigroup generated by — on is bounded and 



analytic, and that it extents consistently to a contractive semigroup on for all q G [1, oo]. 
Now the result is a consequence of |Lam| Corollary. 1.1]. □ 

In order to include lower order terms into the boundary and interface conditions we need 
some preparation. 

Proposition 3.5. Let p G (1, oo) and 9 G (0, 1) be such that d — 1 < 9 p. Then one has 
Dom((^y 9 )c L°°. 



Proof. Since Dom((<^4 p ) e ) = Dom((A p + l) e ) by Theorem 2.2! 
4.1.11], it suffices to show that (A p + iy e maps W into L°°. 
is shown that 



in 



In 



and |Lun[ Lemma 
Pazl Section 2.6] it 



(A p + 1) 



m Jo 



r -1 e -t St dt. 



Now the assertion follows from the estimate of Lemma 2.1S|[r 



□ 



Corollary 3.6. Suppose p G (-^-,00). Then (L p , Dom(^A p )) a _ 1 s continuously embeds 
mtoh 00 . 



Proof. Fix 9 G (— ,« - 1). Then Proposition 3.5 yields Bom((qA p ) e ) C L°°. But 
(L p ,DomM P )) Q .i s C (L p ,Dom(^ p )) eil C [I/, Dom(qA p )] e 

s ' 

by |Lun| Propositions 1.1.3, 1.3.2 and Corollary 2.1.8], and the latter interpolation space 
equals Dom((A p ) 9 ) by Theorem 2.22 iii) □ 



Definition 3.7. Fix b G L P (T U E; dp). Define the operator B : L°° -»■ L p by 

B(f n ,f a ) = (0,bf d ) 
for all / = (/„, /a) e © L p (r U S; dp) S L p . 
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Note that b is allowed to be complex valued. 

Theorem 3.8. Letp G (d— l,oo). Then the operator q(A p + B) satisfies maximal parabolic 
regularity on L p . 

Proof. One deduces from Corollary |3.6 that the operator qB acts continuously on an inter- 
polation space between ~Dom.(qA p ) and L p . Then the result follows from the perturbation 
theorem |Dorl Theorem 6.2]. □ 

Remark 3.9. In a somewhat more general concept B may also depend explicitly on time, 
see jACFP] . 



Now we are in the position to solve the parabolic problem (1.1 )— ( 1.5 ) in terms of the 
realization of the operator A p . 

Theorem 3.10. Let T G (0, oo) and set J = (0, T). Let p G (d - 1, oo), s G (1, oo) and 

a G (-, 1]. Let Q be a bounded domain in IR d with d > 1, letT be an open part of its boundary 
dVL and Ecfl. Adopt the Assumptions 2J2 1 2A and 2.12 Let e G L°° be a positive function 
with a positive essential lower bound and let b as in Definition 3/7, Then the initial value 
problem ( 1.1 ) — ( 1.5 ) admits a solution in the following sense: for all f G L s a (J;h p ) and 
u G (L p , Dom(yLp)) Q _i there is a unique function u G W^' S (J;L, P ) D L s ( J; Dom(A,)) 

s ' 

satisfying 

ev! + A p u + Bu = f, m(0) = u . (3.2) 



Proof. One reformulates (3.2) as 

u + s' 1 A p u + e' 1 Bu = e" 1 f ', u(0)=u . 

Obviously, e~ l f satisfies the same assumptions as /. Moreover, one has Dom(A p ) = 
Dom(e _1 Ap) = Dom(£: -1 (Ap + B)), with equivalent norms. This implies that 

(I/, Bom(A p )) a _ Ks = (L p , DomOT 1 ^ + B))) a _ Ks . 

s 7 s ' 

The assertion then follows from Proposition |3.2| and Theorem |3.8| □ 

Remark 3.11. In the situation of the theorem, the solution depends continuously on 
the data, due to (3.1). Proposition 3.3 gives further regularity properties of a solution. 
Moreover, again by (3.1 ), it is straightforward to see that the solution depends continuously 
on the function e, with respect to the L°°-norm. 

Remark 3.12. Since the coefficient function \i is real valued, the resolvent of qA p commutes 
with complex conjugation on the spaces L p . The latter is also true for the semigroup 
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operators e~ tqAp . Thus, the restriction of <,A p to real spaces also satisfies maximal 
parabolic regularity. If B is induced by a real valued function, then the same is true for 
the operator <,(A P + B). 

Remark 3.13. At the end of this section, let us give more detailed, partly heuristic 
arguments what the real advantage is of the treatment of our parabolic equations in the 
spaces L p . 

When considering the solution u of a parabolic equation u' + Au = f on a Banach space 
X one can form the dual pairing with any if) G X* to obtain 





Of 



(u^) + (Au^) = (f^). 



(3.3) 



E.g., if X = W~ l ' 2 (Q), then one can choose for ip as any element of W ' 2 {fl), but not an 
indicator function of a subset of Q. In our setting, the situation is different: if X = LP, then 
the dual pairing with the indicator function xu of a measurable set U C Q is admissible. 



Then (3.3) reads, there A taken as the L 2 -realization of A 2 , 

u(dx + dp)+ / (A 2 u) (dx + dp) = / f(dx + dp). (3.4) 
Ju Ju 



d_ 

dt 



u 



Since A 2 u G L 2 for almost every time point t we are now at least in principle in the 



position to rewrite j v (A2u) (dx + dp) as a boundary integral and thus to recover from (3.4) 



the 'original' physical balance law for (1.1 )— (1.5) 



Indeed, applying (2.4) with v G C£°(Q), it follows that the distributional divergence 
of pVu is given by the finite Radon measure induced by (A 2 u\q, A 2 u\y) G L 2 {VL) x 
L 2 (£; dl-Ld-i) with respect to dx + dWd-i (see also Remark 2.5). Under certain further 



assumptions on pVu or U one can apply the generalized Gauss-Green theorems of e.g. 
|CTZlj . |Fug| and |Zie] to obtain 



(A 2 u) (dx + dp) 



u 



v ■ pVudT-Ld-i, 



(3.5) 



au 



where v ■ pVu G L 1 (dU ; dHd-i) 
flux', see ibidem. 



is 'the generalized normal component of the corresponding 



Substituting (3.5) in (3.4) gives the desired balance law, as is classical when V ■ pVu is 



an L 2 (f2)-function; compare |Soml Chapter 21] and |ChaL] . As already mentioned in the 
introduction, this is the basis for local flux balances, which are crucial for the foundation 
of Finite Volume methods for the numerical solution of such problems, compare [BRFJ, 
[FL] and [Gar] . 
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4. QUASILINEAR PARABOLIC EQUATIONS 



In this section we treat a nondegenerate quasilinear variant of (1.1 )— ( 1.5 ), including 
nonlinear terms in the dynamic equations on T and £, i.e., 

ed t b{u)-V ■ na{u)Vu = F n (t,u) in Jx(0\S), (4.1) 

u = on J x (0ft (4.2) 

ed t b(u) + v ■ (ia(u)Vu = F r (t, u) on J x T, (4.3) 

sd t b(u) + [vz ■ [ia(u)Vu] = F^(t,u) on J x S, (4.4) 

M (T ) = u in ft U T, (4.5) 

where J = (To, Ti) is a bounded interval. Interesting examples for the nonlinearities on the 
left-hand side are e.g. when b and a are an exponential, or the Fermi-Dirac distribution 
function Jr/2, which is given by 

•Fi/ 2 (s):=4= / i i~4=I d £- 



% Jo 1 + e ? 

Further, in phase separation problems a rigorous formulation as a minimal problem for the 
free energy reveals that a = b' is appropriate. This topic has been thoroughly investigated 
in (§ua|, |QRVl , [GLT] . and [GL2] . see also [GS] and [GrTT] . 

We consider from now on the real part of the spaces L p and the operators A p . For 
simplicity we write L p for L^. As in the linear case we give the quasilinear equation a 
suitable functional analytic formulation, and within this framework the problem will then 



be solved (see Definition 4.2 and Theorem 4.5 below). Again throughout this section we 
fix the numbers 

1 < s < oo and - < a < 1. 

s 



We impose the following conditions on the coefficients on the left-hand side of (4.1)-(4.5). 

Assumption 4.1. The coefficient matrix \i is real- valued, b G W^^°(M) is such that b' is 
positive, and a G W^^°(R) is positive and satisfies J °° a(£) d( = oo = J ^ a(£) d(. 

Note that we do not require monotonicity for a. In particular, terms of the form a(u) = 
7] + \u\ m with 1] > and m > 1 can be treated, that arise e.g. as a regularization of the 
porous medium equation. 

It is in general not to expect that the domain of the realization of —V • fia(v ) V on L p as 
in Section 2.3 is independent of v G L°°(ft). Consider, e.g., the case of a smooth geometry 
with fia(v) equal to a constant on the one hand and a nonsmooth fia(v) on the other hand. 
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This observation motivates our definition of a solution of (4.1)-(4.5), which we describe in 
the following. We put 



-f%(o d(, if e < o. 



Then the assumptions on a imply that 



in Jx (0\E), 
on J x (9fi\r), 
on J x r, 
on Jx S, 

in nur, 



£: R ->■ R is bijective, iT 1 G W 7 "^ 00 (R) , = a, and £(0) = = iT^O). 

In the sequel we identify the functions b, .ft, &T 1 with the Nemytzkii operators they induce. 
The reformulation of (4.1)-(4.5) is based on the so-called Kirchhoff transform w = &(u). 
This (formally) gives a(u)Vu = Vw and d t (b(u)) = ^ (&~ l (w))d t w. Since £(0) = 0, the 
problem (4.1)-(4.5) thus transforms into 

d t w — rj[w) V • //Vw = ri(w)Fn(t, £ _1 (w)) 
w = 

<9 t w + T]{w)v ■ //Vw = rj(w)Fr(t, ^ _1 (w)) 
<9 4 w + r/(w)[z/ s • //Vw] = ?7(w)F 2 (t,^ 1 (w)) 
w(T ) = £(w ) 

where we have set 

t](w) := e' 1 — £ -1 (w). 
For all i G J, let us further define the operator 

VwyFo^JrVln)) onfi\S, 
i?(t,w) := <( ^(w^Fr^^^lr)) on T, (4.6) 
^(w| s )F s (t,^- 1 (w| s )) onE, 

acting on real- valued functions defined on HUT. 

Definition 4.2. Let p G oo), and let A p be the realization of —V ■ //V on L p as in 

Section |23) We say that u G C{[T , 7\]; L°°) is a solution of (fjTHO) on J if 

J*(it) G W l a s {J-U) nL s a (J;Dom(A p )), 

and if w = satisfies 

d t w + r/(w)A p w = R(-,w) on J, w(T ) = &(uq). (4.7) 
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If &(u) is as above, then u G C([T , Ti]; L°°) is already a consequence of Proposi- 
tion |3.3[ Corollary |3.6 and the regularity of Proposition |3.3 shows that in fact u G 
C 7 ([To, Ti\; L°°) for some 7 > 0. For specific choices of .ft additional regularity may carry 
over from R(u) to u. In any case one has u(t, ■) — > uq as t — > To in the L°°-norm. 

Observe further that in the definition it is necessary that &(uo) G (L p , Dom(A p )) a _ 1 . 
It would be interesting to find another description for this condition for a class of non- 
linearities a. If a is constant, then a solution in the above sense can be defined for all 
u G (I/,Dom(i p )) a _i,. 



If = b', then (4.7) is in fact a semilinear problem. This is in particular the case for the 
phase separation problems from above. 



To solve (4.7) we intend to use the following abstract existence and uniqueness result, 
which is proved in |Priij for the temporally unweighted case a = 1. The proof in |Priij 
literally carries over to the weighted case a < 1. 

Proposition 4.3. Let X, D be Banach spaces such that D embeds continuously and densely 
into X. Assume A: (X,D) a _i — > C(D,X) and 71: J x (X, D) 1 — > X are such that 

s ' s ' 

lZ(-,Wo) is measurable for all wo G (X,D) a _i s , that 7Z(-,0) G L s a (J;X) and that for all 

s ' 

M > there are Cm > and r M G L s a (J) with 

\\A(wi) - A(w 2 )\\c(D,x) < C M \w\ - w 2 \\(x,D) x 

and 

\\K(t,Wi) -K(t,w 2 )\\x < r M (t) \\ivt - w 2 \\( X ,d) 1 for a.e. t e J, 

for all wi,w 2 G (X, D) a _i s with 1 < M and \\w 2 \\(x,d) 1 < M. Assume 

further that for any Wq G (X, D) 1 the operator A(w Q ) with domain D on X satisfies 

s 1 

maximal parabolic regularity. 

Then for all wq G (X, D) a _i s there are T* G (To, T\] and a unique maximal solution w 

s ' 

of 

w' + A(w)w = 1Z(-, w) on (T ,T*), w(T ) = w , 
such thatw G W^' S (T , T; X) fl L s a (T , T; D) for all T G (T ,T*). 



4.1 



Let p G (fi,oo 



3.6 



We apply this result to ( 4.7[ ). Suppose b and o satisfy Assumption 
X = L p , D = Dom(A p ) and ^4(w) = rj(w)A p for all w G (L p , Dom(v4^)X Q _i . Corollary 

s ' 

implies that 

(L^Dom(i p )) a _i s cL°°. (4.8) 
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Thus if w G (L p , Dom(A p )) a __ 1 s and ||wo||(LP,Dom( J A } ,)) 1 < M for a given number M, 
then it follows from ( |4.8 ) that the image of f2ur under wq is almost everywhere contained in 
a compact interval that only depends on M. In particular, this gives 7](wo),f](wo)~ 1 G L°°, 
and the operator A(wq) with domain Dom(v4 p ) on L p satisfies maximal parabolic regularity 
by Theorem 3.4| 

The function rj is locally Lipschitz continuous on R. Therefore 

||^(wi) - ^.(w2)||£(Dom(A p ),LP) < \\v( w l) ~ V( W 2) |k°° 

< CW||Wi -W2||l°° < - W 2 ||(LP,Dom(Ap)) _i 

for all Wi,u> 2 G (L p , Dom(A p )) a _ i with ||wj||(LP,Dom(A p )) _i < M for all j G {1,2}. This 
verifies the conditions of the above proposition for A. 

We next present sufficient conditions for the functions Fq, Fr and F-% such that the 



operator R, defined in (4.6), satisfies the conditions for 1Z in Proposition 4.3 



Assumption 4.4. For all £ G R the mappings J — >• R, F r (-,£): J — >• R and 

Fy,{-, £) : J — )■ R are measurable. For all M > there is ru G L s a (J) such that 

|^(t,a)-^n(*,6)| <r M (*)|6-6l 
for a.e. i G J and with |£i|, |£ 2 | < M; and analogous conditions for Fr and Fz. 



Under the above assumption, (4.8) implies that R(-,wq): J — > L p is measurable for 
all wo G (L p , Dom(A,)) a „ i and that -R(-,0) G L s a (J). We verify the Lipschitz prop- 

s 1 

erty for the first component of R. If M > 0, and Wi,W2 G (L p , Dom(A,)) a _ i with 
||wi||(LP,Dom(Ap)) _i < ^ an d 11^2 1| (Lp,Dom(A p )) _i < M, then for a.e. t G J we have 

ll^iln^n^^KIn))-^^)^^^- 1 ^^))!!^) (4.9) 

< \\v( w i\n) - /7( w 2|n)IU->(n)||^h(*,^ -1 (wi|n))||LP(n) 

+ II^I^IU-cnjII-Pn^^KIn)) - •Ph(<,^" 1 («'i|n))||i>(n) 

< Cm (||wi|n - w 2 \n\\L^(n) + rM(0ll# _1 (wi|n) - £~ 1 (w2k)l|.LP(n)) 

< Cm(1 + - W 2 \\(LP,Dom(A p )) _ 1 , 

for a suitable r^/ G L s a (J). The same arguments apply to the other components of R, and 
thus R is as desired to apply the proposition. 
We have proven the main result of this section. 



3.10 



Theorem 4.5. Let p G (^zj,oo), and suppose thattt, T, E, and e are as in Theorem 
that n, b and a are as in Assumption 4.1, and that f, g and h are as in Assumption 
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Then for all u G L°° with ft(u ) G (L p , Dom(A p )) a _i there are T* = T*(u ) G (T ,Ti] 

s ' 

and a unique maximal solution u G C([T , T*); L°°) 0/ (4.1)-(4.5) m £/ie sense of Defini- 
tion 4.2 , 77ms means that for all Tq < T < T* we have 

£(u) G W^ S (T ,T;I/) n (T ,T; Dom(A p )), 

and ^(w) zs i/ie unique solution of 

d t w + r](w)A p w = R(-,w) on(T ,T), w(T ) = £(n ). (4.10) 

Instead of Fq, Fp and F^ one can easily find also non-local maps such that the corre- 
sponding operator R satisfies the condition of Proposition 4J3 One can take for example 
(linear or nonlinear) integral operators with suitable kernel properties. Moreover, in our 
example, Fq maps L°°(Q) into itself, while Fp maps L°°(r) itself, and correspondingly also 
for Fs, i.e., the mapping R has no crossing terms. This is also not necessary in general. 

The nonlinearity in the elliptic operator may also be a nonlocal operator. This case 
arises e.g. in models for the diffusion of bacteria; see [CCj . |ChiLj and references therein. 



We end this section with some comments on the case when (4.1)-(4.5) is semilinear, i.e. 
when b = ^ = id, such that u itself solves the realization (4.10) of the problem. 
The following is a useful criterion for the global existence of solutions. 



Proposition 4.6. Adopt the assumptions of Theorem |4.5[ Suppose in addition that b 
& = id, and let u G C([Tq, T*); L°°) be the maximal solution of (4.1)-(4.5). If 



limsup \\u(t, •) 



(0) 



< OO, 



then u is a global solution, i.e., T* = T\ and u G W^' S (J; L p ) D L s a (J; Dom(A p )) 



Proof. By Proposition 3.2 , for all T < T* the solution u satisfies 

\\u'\\l^(T ,T;Lp) + |MU«(To,TiDam(i4p)) < c(||«o|| (X,Dom(A p )L 1 + \\R(; u) (T ,T;Lp)) , (4.11) 

where c is uniform in T. Observe that \\u(t, Oik 00 < \\ u (t, ') for almost every t by 
the definition of the trace (see Section 2.2). Hence M = ||ii||,L°°(To,T*;L°°) < oo. Estimates 
as in (4.9) yield 

\\R(-,u)\\Li t (To,T";U>) < \\R(-,0)\\li x (T ,T*;U>) + C M (l + \\r M (t ,T*)) • 



Therefore the terms on the left-hand side of (4.11) are bounded uniformly in T. By |Priil 
Corollary 3.2], this implies T* = T v □ 

We finally comment on the asymptotics of solutions. 
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Remark 4.7. Under the assumptions of Theorem 4.5[ in the autonomous semilinear case 
the solutions form a local semiflow in the phase space Dom(v4p), where 6 is sufficiently 
close to 1. Since the resolvent of A p is compact by Lemma 2.1£|[iii) the solution semiflow 
is compact, and bounded orbits are relatively compact. This property is very useful in 
studying the long-time behaviour of solutions. 



5. Concluding remarks 



Remark 5.1. The realization of (1.1 )— ( 1.5 ) in Section |3| still enjoys maximal regularity if 
one adds a term bu in the dynamic equation on J x £ and imposes suitable conditions on 
b. 

Remark 5.2. Everything can be done also for systems which couple in the reactions. 



Remark 5.3. The fundamental result of Priiss (Proposition 4.3) allows to treat the quasi- 



linear problem (4.1 )-(4.5) also in the case where the nonlinearities b and a depend explicitly 



on time. We did not carry out this here for the sake of technical simplicity. 

Remark 5.4. If one requires Q to be a Lipschitz domain and, additionally, imposes a 
certain compatibility condition between T and its complement in the boundary (see [GroJ, 
[HMRSJ), then (-V'/iV + l)" 1 maps Wf ,<? , i.e., the anti-dual space of W r ' 9 , into a Holder 



space, if q > d. If s in Theorem 3.10 /Theorem 4.5 is chosen sufficiently large, then the 



corresponding solutions are even Holder continuous in space and time, compare jDMRTj . 

Remark 5.5. What cannot be treated within this framework is the case where £ moves 
in Q in time. If one wants to include this, the concept of jHaRlj should be adequate, see 
also [HaR3j . 

Remark 5.6. What also cannot be treated within this framework is the case where the 
function e is not away from 0, in particular, if it is on a subset of positive boundary 
measure. This would e.g. affect the case of inhomogeneous Neumann boundary conditions. 
It is known that the resulting problem is of very different functional analytic quality and 
requires different methods, see [Nit] . 
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